function x0 = fn_find_zero(vx, vep)

[a0, a2, a4] = fn_fit_poly(vx, vep);
x0 = sqrt((sqrt(a2^2-4*a0*a4)-a2)/2/a4);

end